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Pair operators for boson and fermion atoms generate 5(7(1, 1) and SU{2) Lie algebras, respec- 
tively. Consequently, the pairing of boson and fermion atoms into diatomic molecules via Feshbach 
resonances, produces 517(1, 1) and 5(7(2) coherent states, making bosonic pairing the matter-wave 
equivalent of parametric coupling and fermion pairing equivalent to the Dicke model of quantum 
, optics. We discuss the properties of atomic states generated in the dissociation of molecular Bose- 

■ Einstein condensates into boson or fermion constituent atoms. The 5(7(2) coherent states pro- 

' duced in dissociation into fermions give Poissonian atom-number distributions, whereas the 5(7(1, 1) 

^SJ I states generated in dissociation into bosons result in super-poissonian distributions, in analogy to 

two-photon squeezed states. In contrast, starting from an atomic gas produces coherent number 
distributions for bosons and super-poissonian distributions for fermions. 
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r— I, I. INTRODUCTION 
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i^H I The behavior of a gas of non-interacting particles close to the absolute zero of temperature depends solely on their 
quantum statistics. Whereas fermions obey Pauli exclusion, manifested in the equal-time anticommutation relations 
of their field operators, bosons are subject to Bose enhancement, implicit in their field operator commutators. For 
' interacting particles, the interaction affects thepair-statistics of fermions and bosons. Pairing models have attracted 
renewed interest since Feshbach resonances Ju Q have been employed to realize molecular Bose-Einstein condensates 
with fermionic [p, and bosonic 0, S, [J constituent atoms, and the ensuing research of the BEC-BCS crossover 

[lOl . [Til . [H, [13II14 lisl. [l6l. Iit\ . Recently, some attention was given to the relation between the quantum statistics 
of the atomic gas in which Feshbach or optical association is performed, and the resulting number-statistics of the 
Q I atomic and molecular fields generated in the process [l^, [l^, [l^, [2l|. It was shown that, whereas the molecular 
field produced in boson association will initially be in a Glauber coherent state, \a), defined such that a\a) — a\a) 
where a is a destruction operator, having constant particle-number fiuctuations [i.e., (An)^/(n) = 1, where n — a^a 
and (An)^ — (n^) — {n)^], the corresponding number distributions for fermions will be super-Poissonian po| with 

■ (An)^ exceeding (n). The coherence of the boson-association field was attributed to collective association, whereas 
l/^ , the chaotic number distributions for fermion-association were related to the individual association of fermionic atom 
00 ■ pairs. 

I Here we explain and quantify these differences in the molecular number-statistics in terms of the commutation 
relations of fermion and boson pair operators. It is well known that pair operators for fermions and bosons generate 
5(7(2) and SU{1,1) algebras, respectively [13, [H, HI]. Consequently, the atomic field produced in the dissociation 
C . of a molecular BEC into fermion atoms will be in an SU (2) coherent state with Poissonian number distribution, 
' whereas boson atoms thus generated, will be in an SU{1, 1) coherent state, corresponding to a squeezed state of the 
J> i Wigner-Weyl algebra, with a super-Poissonian distribution. Using the simple mapping between SU{2) and SU{1, 1) 
it is shown that for association, 5(7(1, 1) coherent states will initially dominate fermion pairing whereas 5(7(2) states 
. ^ will be generated for bosons. Boson association (unlike boson dissociation) is not a collective effect since the molecular 

■ field can be replaced by a macroscopic c-number, rendering the initial molecule production process perfectly linear. 
■ " ' The super-Poissonian statistics of fermion association on the other hand, actually result in from collective behavior 

of the fermionic association. In particular, it will also show up in the degenerate fermionic case, where all atom pairs 
'emit' molecules in-phase. 

In section |ll] we discuss the dynamical equations for associative pairing of fermionic and bosonic atoms. Sec. IIIII 
describes the bosonic and fermionic coherent states of the 5(7(1,1) and 5(7(2) algebras generated by the angular 
momentum like operators respectively. Sec. IIVI discusses the short time dynamics of the dissociation of a molecular 
BEC into bosonic and fermionic atoms, and Sec. fVl concludes the paper. 



II. DYNAMICAL EQUATIONS 



We begin by considering the dynamical equations for associative pairing of fermionic and bosonic atoms, highlighting 
similarities and differences resulting in from the underlying pair statistics of 5(7(2) for fermions and of 5(7(1, 1) for 
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bosons. As shown below, when the atomic motion is slow with respect to other timescales, the atom-molecule 
pairing Hamiltonians map onto two quantum-optical paradigmatic systems; The pairing of fermion atoms is a matter- 
wave equivalent of Dicke superradiance [l^ , whereas the dissociation into bosonic atoms is analogous to parametric 
downconversion [2^. 

A. Fermion atoms - SU{2) algebra 
Consider the single molecular mode association/dissociation Hamiltonian 

^ = H^kc[',Ck,<x+f6^& + .g [6^^Ck,TC_k,x + /i.c.J , (1) 

k,(T V k / 

where ek = h^k^ /2m is the kinetic energy of an atom with mass to, E is the molecular energy, containing kinetic and 
binding contributions, and g is the atom- molecule coupling strength. The annihilation operators for the atoms, Ck.o-, 
obey fermionic anticommutation relations, whereas the molecular annihilation operator h obeys a bosonic commutation 
relation. The single mode approximation is justified when the molecular dispersion due to the presence of a molecular 
momentum spread, is slow with respect to any other timescale in the problem. It becomes exact for a molecular 
BEC, when molecular translation is completely frozen. For simplicity, we have also omitted background non-reactive 
atom-atom scattering. As will be evident from the discussion below, these interactions can be easily incorporated, as 
long as they arc dominated by (k |, — k [) pairing. While this assumption is well-justified for fermions in the BCS 
state, it is a gross oversimplification for bosons. We thus expect that our results will be restricted to the case where 
background open-channel interactions are small with respect to closed-channel atom- molecule coupling, i.e., to narrow 
Feshbach resonances. 

For fermionic atomic field operators, the model Hamiltonian ([1]) can be written using only the atomic SU{2) 
generators [H ^ 

_ „ „ ck _ „t J ck _ 1 / 1 I t ^. . I „t 



'5'_ = Ck,TC^ka , = cLk.|4 , S^= -[-1 + cT_|Ck,T + cLk,^c_ka j , (2) 
obeying the canonical angular-momentum commutation relations 

[5^,^'!] = 25^ [5^5^] = ±5^ (3) 

Using Eqs. Hamiltonian ([T|) may be rewritten as 

i? = ^ ek (2^,^ + 1) + Eb'^b + .9 ^ (fot^k ^ ^k^^ ^ 

k k 

resulting in the Heisenberg equations of motion, 

iS^ = [S^^H] = 2euS^ - 2gS^b, 



ib = [b,H]=£b + g^S^ . (5) 



Defining the Hermitian operators 



equations ([5]) transform into: 

S^^-2euS^^2gByS^, = 2e^S^ - 2gB^S^, = 2g (B^S^ + ByS^) , 

B,=£By-gJ2S^, By^-£B,-gj:^S^ . (7) 
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System ([7|) satisfies the conservation of the individual spin angular momenta, with the SU{2) Casimir operators 

S'^' = S^{S^ - 1) + 5"^ 5*^ = s{s + 1) , (8) 
with s = 1/2, as well as total number conservation 

{Bl + Bl-l/2)+Y,{S^ + m^N/2, (9) 

k 

where N = 2b^b + J2]i,a ^t,a^K<y Defining, 

the dynamical equations ([5]) take the form 

iSX = AuS^ + 2gb^bS^ + gY, S^S^ , 

k' 

iS^ = -Ai,S^-2gb1bS^~gJ2s+S^, 

k' 

iS^ = g{sX- S^) , (11) 

with Ak = £ — 2ek- 

For the degenerate case, ek e, which yields 

Ak ^ A = f - 2e . (12) 

Hamiltonian ([¥]) is, up to an insignificant c-number shift, just the Dicke Hamiltonian [25| and the dynamical equations 
become 

ij+ = AJ+ + 2gb%J^ + gJ+J- , 
ij- = -AJ^ - 2gb'^bJ^ - gJ+J^ , 

iJz = g{J+-J-) , (13) 

where Ji = X)k '^i' ~ Sk'^i'- order to get a closed set of equations for Jj^.J^ and J^, we use the SU{2) 
Casimir operator, 

J,(J,-1) + J+J_ =j(j + l) , (14) 

and number conservation 

b^b^^-ij + J^) , (15) 

where TV = 2b^b + ^j^ ^ cj^ o-Ck.o- is the total number of atoms and j is the number of available energy levels (containing 
at most 4:j particles, because each level can accommodate k t,k j,— k t^^k j atoms). Substituting Eq. (|14p and 
Eq. Uni) into Eqs. we obtain 

ij+ = AJ+-g[3Jf^-iN-2j)J,-f-j-J,] , 
ij- = -AJ^+g[iJl-{N-2j)J,-f-j-J,] , 

ij, = g{J+-J-) . (16) 
Defining the normalized operators 

~ 2(iV/4)3/2 ' ~ 2i(iV/4)3/2 ' ~ iV/4 ' ^ ^ 

we finally obtain the dynamical equations 

- Ajy , (18) 

= -AJ^ + ^ [3J,2 _ (4 _ 2ri) J, - ,72] _ 2|^^^ ^ ^ ^^^^ 



= gVA^J-j, , (20) 
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where 77 = 4:j /N denotes the number of quantum states per particle or the inverse phase space density. For a thermal 
gas T] ^ 1, whereas for a Fermi degenerate gas 77 = 1. For a filled Fermi sea, rj attains its minimal value of unity, and 
Eq. (fT9|l can be replaced by 



jy = -AJ, - ^ [(1 - J.)(l + 3J,)] - ^(1 + J.) . (21) 



B. Boson Atoms - SU{1, 1) algebra 

We next consider the coupling of a molecular BEC into bosonic atom pairs. The single molecular mode Hamiltonian 
reads, 



= ^ ekfl^ k^'T.k + 56^6+17 I ^ ai.ka2,-k + h.c. j 
k,(T \ k / 



(22) 



where ek, £, g, and b have the same meaning as in Eq. ([T]) and the atomic annihilation operators a^r^k, denoting two 
atom species, now obey bosonic commutation relations. 

The pertinent algebra for bosonic atom operators is SU{1,1), because the commutator of = ai.kQ'2.-k and 

Kl = .^aj k is 

[-ft'i , K^l] = (-1 - al,k«i.k - aLka2,-k) = -2K^ , (23) 
so that the three generators K^, K^, obey SU{1, 1) commutation relations 

[K'^,K'l] = -2K^ , [K^, At] = ±K^ , (24) 

differing only in the sign of [A'lf,/ir!!] from the commutation relation between the SU{2) generators, stipulated in 
Eq. ([3]). Hamiltonian ^ is thus replaced by the SU{\, 1) Hamiltonian, 

= ^ fk {2K^ - 1) + Sb^b + gY, {b^K^l + K\b) , (25) 

k k 

leading to the Heisenberg equations of motion. 



iK\ = 




= -2ekKl - 2gb^K^ 


iK^ = 




= 2€kK't + 2gK^b , 


iK^ = 


[K^,H] 


^giKlb-b^K't) , 


ii) = 


[b,H] = 


£b + gY,K^ ■ 



(26) 

For degenerate atomic energy levels, the boson Hamiltonian (j25|) is identical to the model Hamiltonian of parametric 
downconversion (26| . Following the same procedure as in the previous section, we obtain for boson degenerate modes, 

iK.+ = AJC+ - 2gbhK, + gK+K^ , 
ilC- = -AlC- + 2gb^bK^ - gK+K- , 

iK, - .g(/C+ -/C_) , (27) 

where K± — J2k^±^ ^'z ~ Sk-^z' ^+ — I^+b, and JC- = b^K-. In contrast to the unitary SU{2) case where we 
had 

-3 < {Jz) < -J+mm{N/2,2j}, 
reducing to ~j < (J^) < j for N = 4j, we now have 

k < {Kz) <k + N/2, 
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where 4fc denotes the number of boson atomic modes. However, we can still eliminate K± and h by using number 
conservation and the SU{\, 1) Casimir operator, 

6t6 = y - {K, - k) , (28) 
K,{K, - 1) - K+K- = k{k - 1) , (29) 

resulting in the dynamical equations 

ilC+ = A/C+ + g [iKl ~{N + 2k)K^ ^ + k - K,] , 
i)C- = -A/C_ - g [S/if -{N + 2k)K, -k^ + k~K,] , 

ik, = .g(/C+ -/C_) . (30) 

We define, as we did for fermion atoms, 

_ /C+ +/C- _ K.+ -K.- _ 

" ~ 2(Ar/4)3/2 ' '^^ ~ 2i(iV/4)3/2 ' ~ (7V/4) ' ^ 

and using these definitions, the dynamical equations are transformed to the final form 





- A/C, , 






= A/C,. {iKl (4 + 2^)/C, ,7^]- 











(32) 

In order to gain better insight on the relation between the fermion equation (PO)) and the boson equation ([5^ we 
define the number difference operator Lz = (27ib — Ua) /N = 1 + — /C^, whose expectation value, like the expectation 
value of J7^, corresponds to the atom- molecule population imbalance. With this definition, we have 

ilCl - (4 + 27?) JCz-v^^ 3Cl - (2 + 47/) - 1 , (33) 

/C, - 7/ = 1 - , (34) 
and the dynamical equations (|32p assume the form 

4 = -A/C, - ^ (3£, + 1) (£, - 1) + ^ [(4fc - 1)£. + 1] , 

^ V -'V 



£2 = -gVNICy . (35) 
The two atomic-modes case with Hamiltonian 

H = e (^a\ai + 4a2^ + Sb^b + g (b^aia2 + 44^) (36) 

is obtained from Eqs. (|35p by substituting k — 1/2 (because the minimum value of K^, obtained where no atoms are 
present, is 1/2). It is easily verified that the resulting equations of motion for ICx,ICy,Cz are identical up to the sign of 
with the fermion equations for Jx^Jy-, Jz when 77 = 1. Noting that Jz for 77 = 1 and Lz have inverse interpretation 
(i.e. the former equals {jia — 2nb)IN and the latter is (277,;, — na)/N) we see that the dynamics of degenerate fermion 
association maps into two-mode boson dissociation and vice versa. 



III. COHERENT STATES 



Having developed the time-dependent many-body formalism and established the connection with the quantum- 
optical paradigms, we turn to the investigation of the dissociation of a molecular BEC consisting either of fermionic 
or bosonic constituent atoms. For sufRciently short times, we neglect molecular fluctuations and treat the molecular 
field b as an undepleted pump, replacing it by the c-number N 12. The resulting Hamiltonian for fermion (boson) 
atoms under this approximation, thus consists of linear sums of operators generating the SU{2) {SU{1, 1)) algebra. 
Consequently, generalized coherent matter states of the pertinent Lie algebras [2^ . can be dynamically generated in 
the dissociation of molecular BECs. In this section we briefly discuss the properties of SU{2) {SU{1,\)) coherent 
states generated in the dissociation of a molecular BEC into fermion (boson) atoms. 
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A. SU{2) Coherent states 

The generalized coherent states associated with the unitary representations of the SU{2) Lie algebra, are 
parametrized by the two polar (Euler) angles 6 and corresponding to rotations of the fully stretched atomic vacuum 
state |j, —j) (where | j, m) denote the usual mutual eigenstates of the Casimir operator and of the number differ- 
ence operator J^, i.e., J^|j,?7i) ~ j{j + m), J^jj, m) = m\j,m) with m — — j, . . . ,j) about the and Jz axes, 
respectively: 



16*, 0) = cxp {-i(j)Jz) exp {-i6Jx) \j, ~j) = cxp {aJ+ - a*J^) |j, -j) 



(37) 



with a ~ (6*72) exp(— i(^) [23|. Definition (|57|) results in the familiar expansion of SU{2) coherent states in terms of 
number (Fock) states 



tan" - 
2 



E 



tan ( - ) exp(-^ 



i+m 



2j 

i + TO 



1/2 



(38) 



Using either Eq. pT]) or Eq. pS]) it is easily verified that 



(Jy) 



J smf/ cos ( 
J sin 9 sin (j 
j cos 6*, 



(39) 
(40) 
(41) 



so that the expectation values of J are restricted to the Bloch sphere of radius j, as depicted in Fig. [TJ The coherent 
state variance of these operators is 



AJ2 



I (1 - sin^ d cos^ i 

^ (1 - sin^ 9 sin^ ,; 

^ sin^ . 
2 



(42) 
(43) 
(44) 



The total variance of coherent states is thus also bound because |AJp = 
commutation relations lead to the uncertainty relations 



(J'>-(J)' =jU + l)-f = 3- The 



1, 



(45) 



where 



eji are the SU (2) structure constants. In particular, for and Jy we have 



AJ,AJ^> i|(J,)|, 



(46) 



In Fig. [T] we plot the expectation values of J/j = {u,v,w) for SU{2) coherent states, as well as the AJ^ and AJy 
variance of ten such states. Coherent states for which inequality (|46p is an equality are referred to as 'intelligent 
states' or 'ideal coherent states'. From Eqs. and (|^T|) we obtain that SU{2) intelligent states are found 

for = 0, 7r/2, tt, 37r/2 and arbitrary 0, as depicted by dashed curves in Fig. [TJ A subset of the intelligent states are 
the minimum uncertainty states with ~ 0, 7r/2, tt, 37r/2 and = tt/2 (denoted by magenta ellipsoids in Fig. [1]), for 
which the r.h.s. of Eq. (|46p is minimized, with AJ^AJy = 0. While the states with 6 — and arbitrary (yellow 
disks) are also intelligent, their value of AJ^AJy = j/2 is in fact maximal and larger than AJ^AJy = j/4 obtained 
for the non-intelligent states denoted by cyan disks. 



B. SU(1,1) Coherent states 



The mutual eigenstates of the SU{1, 1) Casimir operator ((29)) and of Kz form the basis set: 



(47) 
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FIG. 1: Bloch sphere (shaded sheU) and coherent states of SU{2). Dashed black curves mark inteUigent coherent states. EUipses 
denote AJx and AJy variance for ten coherent states: the atomic and molecular vacuum states (yellow), four non-intelligent 
states (cyan) and the four squeezed, minimum uncertainty states (magenta). 



Kz\k, n) = {k + n)\k, n), 
with n = 0, 1, . . . , N/2. In analogy to Eq. ((37|) . SU{1, 1) coherent states are obtained as 

=exp(/3A%-/r/v_)|fc,0) , 



(48) 



(49) 



with [3 = —{9/2) cxp(— [l^. Power-series expansion of the exponents in Eq. ([49]) gives the SU{1, 1) coherent states 
in terms of the number states |fc,n), 



\0A) - 



1 - tanh^ 



k N/2 

E 

)i=0 



tanh - ) exp(— i 



T{n + 2k) 
n\T{2k) 



1/2 



Consequently, the expectation values of K are 



{Kx) ~ A:sinh6'cos( 
{Ky) — kiimhO sa\(f 
{K^ = /c cosh 61, 



(50) 



(51) 
(52) 
(53) 



so that the motion of the vector (K) is restricted to the paraboloid {Kz)"^ — {K^)'^ — (Ky)'^ = fc^ (see Figl2]), as could 
be expected from the SU{1, 1) Casimir in Eq. (|29p . The variance of the SU{1, 1) generators for the coherent states 
([50)1 are given by 



A'f 

AKl = — (1 + sinh^6lcos2, 



M 

~2 

M 



1 + sinh^ 6'sin^ ( 



■ sinh^ d . 



(54) 
(55) 
(56) 
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FIG. 2: Surface of motion of (K) — {u,v,w) (shaded paraboloid) and coherent states of SU{1, 1). Dashed black curves mark 
intelligent coherent states. Ellipses denote the AKx and Aivy variances for nine intelligent coherent states. Whereas the atomic 
vacuum (yellow ellipse) is a minimum uncertainty state with equal variance in the ICx and ICy directions, other intelligent states 
(magenta,blue) are squeezed. 

Due to the possibility of multiple occupation in any single mode, neither the expectation values nor the variance of 
the K operators are bound. Since the structure constants of the two algebras differ only in sign, the uncertainty 
relations of SU{1, 1) are the same as for SU{2), e.g., 

AK^AKy>^\{K,)l (57) 

and we can define intelligent and minimum-uncertainty states as we did for SU(2) in the previous subsection. In Fig 
[5] we plot the expectation values of K for SU{2) coherent states, as well as the AK^ and AKy variance of nine such 
states. It is clear from Eqs. ([M)) . ([55]) . and ((55]) that the intelligent states will be obtained for = 0, tt/2, tt, 3tt/2 and 
arbitrary 0. 

C. Generalized Squeezing 

As is clear from Eqs. (^5)1 and ([57)) . the minimum fluctuation product of these two observables, depends on the 
expectation value of the remaining generator: 

AXf<^\c^^{X,)\, and > 1|4(X,)|. (58) 

Generalized squeezed states of Lie algebras generated by o ..^ are defined as those states for which the variance 

in one observable has been reduced at the expense of another . 

It is clear from Eqs. (|39p - (|44p that starting from a fermionic atomic vacuum and inducing a rotation about the 
Jx axis (e.g., by choosing A = and g = g*), the variance in Jy will be squeezed at the expense of Jx fluctuations 
(see Fig. [T]) because 1 — sin^ 6* < |cos0| for all 6. Similarly, rotating about the Jy axis (e.g., by choosing A = 
and g = —g*) will result in generalized squeezing of AJ^- The same is also true for rotations of a boson atomic 
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FIG. 3: (color online) Number distributions as a function of rescaled time in the dissociation of 500 molecules into fermion (a) 
and boson (b) constituent atoms. 

vacuum, as seen from Eqs. (|5ip - (|56p . However, the way to attain generalized squeezing for fermions and bosons is 
quite different. Whereas with fermion atoms, squeezing in is obtained by reduction of its fluctuations, keeping a 
fixed Jy variance (thereby reducing the product AJa;AJy = IJ2I/2), the same goal is attained for boson atoms by 
increasing the variance of Ky (thereby increasing AK^AKy = |i4r2|/2)and keeping fluctuations fixed. 

IV. SHORT TIME DYNAMICS 

Since the atomic vacuum state (yellow disk in Figs. [T] and [2]) is a coherent state, the atomic states produced in the 
dissociation of a molecular BEC will initially be coherent states of SU{1, 1) if the constituent atoms are bosons or 
of SU (2) when the constituent atoms are fermions. There should thus be a significant difference in the short-time 
dynamics and in the initial fluctuations between the two cases. For bosons, one expects exponential amplification 
of atom-number and atom number fluctuations, whereas fermion number growth is more moderate and fluctuations 
remain bound. Physically, the source of these differences is in the underlying mechanisms of Bose-stimulation of 
dissociating boson pairs, leading to the dynamical instability of the atomic vacuum, and Pauli blocking of dissociating 
fermion pairs. 

In order to verify the formation of such coherent states and generalized squeezing, we have carried out many- 
particle simulations of molecular BEC dissociation into either fermionic or bosonic constituent atoms. For sufficiently 
short propagation times the molecular field is to a good approximation undepleted, and the generalized operators 
J'ijICi coincide, up to an insignificant c-number, with the SU generators Ji,Ki. The atomic states during the initial 
stage of dissociation are thus approximately SU{2) and SU{1, 1) coherent states respectively for fermion and boson 
constituents. In what follows, we shall numerically investigate to what extent do the generalized coherent states and 
squeezed fluctuations depicted in section IIIII for an undepleted pump approximation, carry through to the operators 
J7i,/Ci, which account for pump depletion and fluctuations. In Fig. [31 we plot the atom-number distribution as a 
function of the rescaled time r = g^/Nt. The dissociation into fermion constituents shown in Fig. [3^, exhibits a 
Poissonian atom- number distribution with bound fluctuations as expected from Eqs. (|4ip and (|44p . Dissociation 
into boson constituents however, results in a super-Poissonian number distribution with an exponential growth of 
fluctuations as predicted in ([55)1 and ([55)) . 

Generalized SU (2) squeezing and its extension into the depleted-pump regime is illustrated in Fig. [4| and Fig. [5| 
where the variance in J^x and Sfy in the dynamical evolution of the (fermion) atomic vacuum state, are plotted as a 
function of time. In Fig. [31 the phase of the association pump is 1^9 = 0, corresponding to rotation about the u axis of 
the Bloch sphere of Fig. [1] i.e., = 0. The coalescence of the variance product AJ^AJy with the expectation value 
of the commutator {[Jx^Jy]) demonstrates that the generated SU{2) coherent state is indeed intelligent. Squeezing 
of the AJy variance while keeping a fixed ASx is observed, in agreement with the undepleted pump prediction. This 
reduction of jTj, fluctuations eventually results in the expected minimum uncertainty state with AJ^AJy k, 0. In 
comparison, the evolution of variances for a (/? = 7r/4 phase of the pump (corresponding to rotation along the = 7r/4 
circle on the Bloch sphere of Fig. [T|) is shown in Fig. [51 The SU{2) coherent states produced during this evolution, 
are non- intelligent, with equal AjTx and AjTj, variances whose product is larger than the uncertainty limit. 

The time evolution of variances in the propagation of the atomic vacuum state for boson atoms, with (p = 0, 7r/4, 
is shown respectively in Fig. [6] and Fig. [71 Unlike the fermion case where fluctuations are bound, we observed 
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FIG. 4: (color online) Quadrature variances AJ'^ (dashed blue line), AJ'y (solid red line), and AJ^AJy (dash-dotted green 
line) as a function of the rescaled time g\fNt, in the dissociation of a molecular BEG made of fermion-dimers. Pump phase 
is = 0. The solid magenta line denotes the uncertainty limit \i\j2:,Jy\)\l1 which can be approximated at early times to be 
2ibpiJz|/(iV/2)2 (dotted black curve). 




resettled time 

FIG. 5: (color online) Same as Fig. |3l with pump phase = 7r/4. Quadrature variances AJ^, ^^v^ AJ^AJy are all equal 
and greater than the uncertainty limit |([^7m7y])|/2 or its short-time approximation 2|&j^|^2|/(Af/2) (dotted black curve). 

rapid increase in the A/Cj, variance with a fixed A/C^,, corresponding to motion on the = parabola in Fig. [21 
While the variance product A/C^^A/Cy grows exponentially with time, its initial evolution traces the uncertainty limit 
[([/C^^, /Cy]) /2 indicating that the produced states are indeed 5'f7(l, 1) intelligent coherent states. Here too, propagation 
with a zero phase of the pump leads to the expected generalized squeezing. For (/? = 7r/4 however, there is no squeezing 
as both A/Cx and A/Cy fluctuations are equal and exponentially growing. The SU(\^ 1) coherent states produced are 
non-intelligent because the variance product A/C^jA/Cj, is larger than the uncertainty limit. 

The agreement between the numerically-exact variance dynamics of Figs. |4][7] and the undepleted pump pictures 
of Fig. [1] and Fig. [21 demonstrates that the different collective dynamics predicted for fermion and boson constituent 
atoms can indeed be interpreted in terms of fluctuations of dynamical variables quadratic in the atomic creation or 
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FIG. 6: (color online) Quadrature variances AJC^ (dashed blue curve), AJCy (solid red curve), and AJCxAJCy (dash-dotted green 
curve) as a function of the rescaled time g\/Nt, in the dissociation of a molecular BEG made of boson-dimers. Pump phase 
is = 0. The solid magenta curve denotes the uncertainty limit \{[K.x,K.y\)\/2 which coincides with 4|6p |/Cz |/(Ai'/2)^ (dotted 
black curve) at early times where molecular depletion is negligible. 
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FIG. 7: (color online) Same as Fig. U] with pump phase — 7r/4. Quadrature variances A/C^, A/C^, AK,xAK,y are all equal 
and greater than the uncertainty limit |{[/Ci,, /Cj,]) |/2 or its short-time approximation 2|&p|A^2 |/(7V/2)^ (dotted black curve). 



annihilation operators. The reduction of fluctuations of these variables is related to the formation of coherent states 
of the SU{2) and ^^^(l, 1) Lie algebras. The same qualitative picture seems to apply to the depleted pump regime. 

Finally, the mapping between fermion and boson dynamics, manifested in the equivalence of Eqs. ()18p - (j20p with 
r] = 1 and Eqs. ()35|1 with k = 1/2, is illustrated in Fig.[8]whcrc atom number distributions are plotted as a function 
of time throughout the association of a Fermi (Fig. [5^) and Bose (Fig. ^p) atomic quantum gas. Fermion association 
is mapped onto boson dissociation (Fig. [8)3) while boson association coincides with fermion dissociation (Fig. [3^), 
generating similar coherent states. 
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FIG. 8: (color online) Number distributions as a function of rescaled time for the association of quantum degenerate gases of 
fermions (a) and bosons (b) with N = 1000. 



V. SUMMARY AND CONCLUSIONS 

We have established a connection between the collective behavior in boson and fermion pairing via Feshbach 
resonances and the generation of coherent states of the SU{1, 1) and SU{2) Lie algebras. This relation provides a 
new viewpoint on the quantum statistics of atom-molecule quantum gas processes. The equivalence of molecular 
BEC dissociation with the Dickc model for fermion atoms p5| and with parametric downconversion for boson atoms 
[2^, known for some time in the quantum optics literature, is put into context as these two quantum systems 
are paradigmatic examples of the aforementioned algebras. The well known squeezing of fluctuations in dynamical 
variables linear in the atomic creation or annihilation operators (e.g. quadrature squeezing) during the dissociation 
into boson constituents, may be viewed as generalized SU{1, 1) squeezing of pair fluctuations, quadratic in the 
atomic creation and annihilation operators. Similarly, the coherent evolution during dissociation into fermion atoms 
corresponds to the generation of a minimum uncertainty SU{2) coherent state. Our numerical simulations indicate 
that the same qualitative picture applies to the fluctuation of operators that account for molecular pump depiction. 
The presentation of the atom-molccule system in terms of these generators offers a link between the fermion- and 
boson-constituent atom cases due to the close relation and direct mapping between the underlying Lie algebras. 
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